Adaptive Selection of Cis-regulatory Elements in the Han Chinese.

Cis-regulatory elements have an important role in human adaptation to the living environment. However, the lag in population genomic cohort studies and epigenomic studies, hinders the research in the adaptive analysis of cis-regulatory elements in human populations. In this study, we collected 4,013 unrelated individuals and performed a comprehensive analysis of adaptive selection of genome-wide cis-regulatory elements in the Han Chinese. In total, 12.34% of genomic regions are under the influence of adaptive selection, where 1.00% of enhancers and 2.06% of promoters are under positive selection, and 0.06% of enhancers and 0.02% of promoters are under balancing selection. Gene ontology enrichment analysis of these cis-regulatory elements under adaptive selection reveals that many positive selections in the Han Chinese occur in pathways involved in cell-cell adhesion processes, and many balancing selections are related to immune processes. Two classes of adaptive cis-regulatory elements related to cell adhesion were in-depth analyzed, one is the adaptive enhancers derived from neanderthal introgression, leads to lower hyaluronidase level in skin, and brings better performance on UV-radiation resistance to the Han Chinese. Another one is the cis-regulatory elements regulating wound healing, and the results suggest the positive selection inhibits coagulation and promotes angiogenesis and wound healing in the Han Chinese. Finally, we found that many pathogenic alleles, such as risky alleles of type 2 diabetes or schizophrenia, remain in the population due to the hitchhiking effect of positive selections. Our findings will help deepen our understanding of the adaptive evolution of genome regulation in the Han Chinese.

Recent positive selection signatures reveal phenotypic evolution in the Han Chinese population.

Characterizing natural selection signatures and relationships with phenotype spectra is important for understanding human evolution and both biological and pathological mechanisms. Here, we identified 24 genetic loci under recent selection by analyzing rare singletons in 3946 high-depth whole-genome sequencing data of Han Chinese. The loci include immune-related gene regions (MHC cluster, IGH cluster, STING1, and PSG), alcohol metabolism-related gene regions (ADH1B, ALDH2, and ALDH3B2), and the olfactory perception gene OR4C16, in which the MHC cluster, ADH1B, and ALDH2 were also identified by TOPMed and WestLake Biobank. Among the signals, the IGH cluster is particularly interesting, in which the favored allele of variant 14_105737776_C_T (rs117518546, IgG1-G396R) promotes immune response, but also increases the risk of an autoimmune disease systemic lupus erythematosus (SLE). It is also surprising that our newly discovered ALDH3B2 evolved in the opposite direction to ALDH2 for alcohol metabolism. Besides monogenic traits, we found that multiple complex traits experienced polygenic adaptation. Particularly, multi-methods consistently revealed that lower blood pressure was favored in natural selection. Finally, we built a database named RePoS (recent positive selection, http://bigdata.ibp.ac.cn/RePoS/) to integrate and display multi-population selection signals. Our study extended our understanding of natural evolution and phenotype adaptation in Han Chinese as well as other populations.

Characterization of genome-wide STR variation in 6487 human genomes.

Short tandem repeats (STRs) are abundant and highly mutagenic in the human genome. Many STR loci have been associated with a range of human genetic disorders. However, most population-scale studies on STR variation in humans have focused on European ancestry cohorts or are limited by sequencing depth. Here, we depicted a comprehensive map of 366,013 polymorphic STRs (pSTRs) constructed from 6487 deeply sequenced genomes, comprising 3983 Chinese samples (~31.5x, NyuWa) and 2504 samples from the 1000 Genomes Project (~33.3x, 1KGP). We found that STR mutations were affected by motif length, chromosome context and epigenetic features. We identified 3273 and 1117 pSTRs whose repeat numbers were associated with gene expression and 3'UTR alternative polyadenylation, respectively. We also implemented population analysis, investigated population differentiated signatures, and genotyped 60 known disease-causing STRs. Overall, this study further extends the scale of STR variation in humans and propels our understanding of the semantics of STRs.

NPInter v5.0: ncRNA interaction database in a new era.

Noncoding RNAs (ncRNAs) play key regulatory roles in biological processes by interacting with other biomolecules. With the development of high-throughput sequencing and experimental technologies, extensive ncRNA interactions have been accumulated. Therefore, we updated the NPInter database to a fifth version to document these interactions. ncRNA interaction entries were doubled from 1 100 618 to 2 596 695 by manual literature mining and high-throughput data processing. We integrated global RNA-DNA interactions from iMARGI, ChAR-seq and GRID-seq, greatly expanding the number of RNA-DNA interactions (from 888 915 to 8 329 382). In addition, we collected different types of RNA interaction between SARS-CoV-2 virus and its host from recently published studies. Long noncoding RNA (lncRNA) expression specificity in different cell types from tumor single cell RNA-seq (scRNA-seq) data were also integrated to provide a cell-type level view of interactions. A new module named RBP was built to display the interactions of RNA-binding proteins with annotations of localization, binding domains and functions. In conclusion, NPInter v5.0 (http://bigdata.ibp.ac.cn/npinter5/) provides informative and valuable ncRNA interactions for biological researchers.

Characterizing mobile element insertions in 5675 genomes.

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

NyuWa Genome resource: A deep whole-genome sequencing-based variation profile and reference panel for the Chinese population.

The lack of haplotype reference panels and whole-genome sequencing resources specific to the Chinese population has greatly hindered genetic studies in the world's largest population. Here, we present the NyuWa genome resource, based on deep (26.2×) sequencing of 2,999 Chinese individuals, and construct a NyuWa reference panel of 5,804 haplotypes and 19.3 million variants, which is a high-quality publicly available Chinese population-specific reference panel with thousands of samples. Compared with other panels, the NyuWa reference panel reduces the Han Chinese imputation error rate by a margin ranging from 30% to 51%. Population structure and imputation simulation tests support the applicability of one integrated reference panel for northern and southern Chinese. In addition, a total of 22,504 loss-of-function variants in coding and noncoding genes are identified, including 11,493 novel variants. These results highlight the value of the NyuWa genome resource in facilitating genetic research in Chinese and Asian populations.

SmProt: A Reliable Repository with Comprehensive Annotation of Small Proteins Identified from Ribosome Profiling.

Small proteins specifically refer to proteins consisting of less than 100 amino acids translated from small open reading frames (sORFs), which were usually missed in previous genome annotation. The significance of small proteins has been revealed in current years, along with the discovery of their diverse functions. However, systematic annotation of small proteins is still insufficient. SmProt was specially developed to provide valuable information on small proteins for scientific community. Here we present the update of SmProt, which emphasizes reliability of translated sORFs, genetic variants in translated sORFs, disease-specific sORF translation events or sequences, and remarkably increased data volume. More components such as non-ATG translation initiation, function, and new sources are also included. SmProt incorporated 638,958 unique small proteins curated from 3,165,229 primary records, which were computationally predicted from 419 ribosome profiling (Ribo-seq) datasets or collected from literature and other sources from 370 cell lines or tissues in 8 species (Homo sapiens, Mus musculus, Rattus norvegicus, Drosophila melanogaster, Danio rerio, Saccharomyces cerevisiae, Caenorhabditis elegans, and Escherichia coli). In addition, small protein families identified from human microbiomes were also collected. All datasets in SmProt are free to access, and available for browse, search, and bulk downloads at http://bigdata.ibp.ac.cn/SmProt/.

Dynamic-BM: multispecies Dynamic BodyMap database from temporal RNA-seq data.

Biological processes, especially developmental processes, are often dynamic. Previous BodyMap projects for human and mouse have provided researchers with portals to tissue-specific gene expression, but these efforts have not included dynamic gene expression patterns. Over the past few years, substantial progress in our understanding of the molecular mechanisms of protein-coding and long noncoding RNA (lncRNA) genes in development processes has been achieved through numerous time series RNA sequencing (RNA-seq) studies. However, none of the existing databases focuses on these time series data, thus rendering the exploration of dynamic gene expression patterns inconvenient. Here, we present Dynamic BodyMap (Dynamic-BM), a database for temporal gene expression profiles, obtained from 2203 time series of RNA-seq samples, covering >25 tissues from five species. Dynamic-BM has a user-friendly Web interface designed for browsing and searching the dynamic expression pattern of genes from different sources. It is an open resource for efficient data exploration, providing dynamic expression profiles of both protein-coding genes and lncRNAs to facilitate the generation of new hypotheses in developmental biology research. Additionally, Dynamic-BM includes a literature-based knowledgebase for lncRNAs associated with tissue development and a list of manually selected lncRNA candidates that may be involved in tissue development. Dynamic-BM is available at http://bioinfo.ibp.ac.cn/Dynamic-BM.

The long noncoding RNA ASNR regulates degradation of Bcl-2 mRNA through its interaction with AUF1.

The identification and characterization of long non-coding RNAs (lncRNAs) in diverse biological processes has recently developed rapidly. The large amounts of non-coding RNAs scale consistent with developmental complexity in eukaryotes, indicating that most of these transcripts may have functions in the regulation of biological processes and disorder in the organisms. In particular, Understanding of the overall biological significance of lncRNAs in cancers still remains limited. Here, we found a nuclear-retained lncRNA, termed Lnc_ASNR (apoptosis suppressing-noncoding RNA), which serves as a repressor of apoptosis. Lnc_ASNR was discovered in a set of microarray data derived from four kinds of tumor and adjacent normal tissue samples, and displayed significant up-regulation in the tumor tissues. Using an RNA-pull down assay, we found that Lnc_ASNR interacted with the protein ARE/poly (U)-binding/degradation factor 1(AUF1), which is reported to promote rapid degradation of the Bcl-2 mRNA, an inhibitor of apoptosis. Lnc_ASNR binds to AUFI in nucleus, decreasing the cytoplasmic proportion of AUF1 which targets the B-cell lymphoma-2 (Bcl-2) mRNA. Taken together, the overall effect of Lnc_ASNR expression is thus a decrease in cell apoptosis indicating that Lnc_ASNR may play a vital role in tumorigenesis and carcinogenesis.

BioCircos.js: an interactive Circos JavaScript library for biological data visualization on web applications.

UNLABELLED: We here present BioCircos.js, an interactive and lightweight JavaScript library especially for biological data interactive visualization. BioCircos.js facilitates the development of web-based applications for circular visualization of various biological data, such as genomic features, genetic variations, gene expression and biomolecular interactions. AVAILABILITY AND IMPLEMENTATION: BioCircos.js and its manual are freely available online at http://bioinfo.ibp.ac.cn/biocircos/ CONTACT: rschen@ibp.ac.cn SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Long Noncoding RNA ADINR Regulates Adipogenesis by Transcriptionally Activating C/EBPα.

C/EBPα is a critical transcriptional regulator of adipogenesis. How C/EBPα transcription is itself regulated is poorly understood, however, and remains a key question that needs to be addressed for a complete understanding of adipogenic development. Here, we identify a lncRNA, ADINR (adipogenic differentiation induced noncoding RNA), transcribed from a position ∼450 bp upstream of the C/EBPα gene, that orchestrates C/EBPα transcription in vivo. Depletion of ADINR leads to a severe adipogenic defect that is rescued by overexpression of C/EBPα. Moreover, we reveal that ADINR RNA specifically binds to PA1 and recruits MLL3/4 histone methyl-transferase complexes so as to increase H3K4me3 and decrease H3K27me3 histone modification in the C/EBPα locus during adipogenesis. These results show that ADINR plays important roles in regulating the differentiation of human mesenchymal stem cells into adipocytes by modulating C/EBPα in cis.

Transcriptome profiling of esophageal squamous cell carcinoma reveals a long noncoding RNA acting as a tumor suppressor.

Esophageal Squamous Cell Carcinoma (ESCC) is among the most common malignant cancers worldwide. In the past, extensive efforts have been made to characterize the involvement of protein-coding genes in ESCC tumorigenesis but few for long noncoding RNAs (lncRNAs). To investigate the transcriptome profile and functional relevance of lncRNAs, we performed an integrative analysis of a customized combined lncRNA-mRNA microarray and RNA-seq data on ESCCs and matched normal tissues. We identified numerous lncRNAs that were differentially expressed between the normal and tumor tissues, termed "ESCC-associated lncRNAs (ESCALs)", of which, the majority displayed restricted expression pattern. Also, a subset of ESCALs appeared to be associated with ESCC patient survival. Gene set enrichment analysis (GSEA) further suggested that over half of the ESCALs were positively-or negatively-associated with metastasis. Among these, we identified a novel nuclear-retained lncRNA, named Epist, which is generally highly expressed in esophagus, and which is down-regulated during ESCC progression. Epist over-expression and knockdown studies further suggest that Epist inhibits the metastasis, acting as a tumor suppressor in ESCC. Collectively, our analysis of the ESCC transcriptome identified the potential tumor suppressing lncRNA Epist, and provided a foundation for future efforts to identify functional lncRNAs for cancerous therapeutic targeting.

Functional Characterization of Long Noncoding RNA Lnc_bc060912 in Human Lung Carcinoma Cells.

Long noncoding RNAs (lncRNAs) are pervasively transcribed in the human genome. Recent studies suggest that the involvement of lncRNAs in human diseases could be far more prevalent than previously appreciated. Here we have identified a lncRNA termed Lnc_bc060912 whose expression is increased in human lung and other tumors. Lnc_bc060912 is 1.2 kb in length and is composed of two exons. The expression of Lnc_bc060912 was repressed by p53. Lnc_bc060912 suppressed cell apoptosis. Using a recently developed method for RNA-pulldown with formaldehyde cross-linking, we found that Lnc_bc060912 interacted with the two DNA damage repair proteins PARP1 and NPM1. Together, these results suggest that Lnc_bc060912, via PARP1 and NPM1, affects cell apoptosis and may play important roles in tumorigenesis and cancer progression.

Distinct MicroRNA Subcellular Size and Expression Patterns in Human Cancer Cells.

Introduction. Small noncoding RNAs have important regulatory functions in different cell pathways. It is believed that most of them mainly play role in gene post-transcriptional regulation in the cytoplasm. Recent evidence suggests miRNA and siRNA activity in the nucleus. Here, we show distinct genome-wide sub-cellular localization distribution profiles of small noncoding RNAs in human breast cancer cells. Methods. We separated breast cancer cell nuclei from cytoplasm, and identified small RNA sequences using a high-throughput sequencing platform. To determine the relationship between miRNA sub-cellular distribution and cancer progression, we used microarray analysis to examine the miRNA expression levels in nucleus and cytoplasm of three human cell lines, one normal breast cell line and two breast cancer cell lines. Logistic regression and SVM were used for further analysis. Results. The sub-cellular distribution of small noncoding RNAs shows that numerous miRNAs and their isoforms (isomiR) not only locate to the cytoplasm but also appeare in the nucleus. Subsequent microarray analyses indicated that the miRNA nuclear-cytoplasmic-ratio is a significant characteristic of different cancer cell lines. Conclusions. Our results indicate that the sub-cellular distribution is important for miRNA function, and that the characterization of the small RNAs sub-cellular localizome may contribute to cancer research and diagnosis.

A differential sequencing-based analysis of the C. elegans noncoding transcriptome.

Noncoding RNAs are increasingly being recognized as important players in eukaryote biology. However, despite major efforts in mapping the Caenorhabditis elegans transcriptome over the last couple of years, nonpolyadenylated and intermediate-size noncoding RNAs (is-ncRNAs) are still incompletely explored. We have combined an enzymatic approach with full-length RNA-Seq of is-ncRNAs in C. elegans. A total of 473 novel is-ncRNAs has been identified, of which a substantial fraction was associated with transcription factor binding sites and developmentally regulated expression patterns. Analysis of sequence and secondary structure permitted classification of more than 200 is-ncRNAs into several known RNA classes, while another 33 is-ncRNAs were identified as belonging to two previously uncharacterized groups of is-ncRNAs. Three of the unclassified is-ncRNAs contain the 5' Alu domain common to SRP RNAs and specifically bound with the SRP9/14 heterodimer in vitro. One of these (inc394) showed 65% sequence identity with the human, neuron-specific BC200 RNA. Structure-based clustering analysis and in vitro binding experiments supported the notion that the nematode stem-bulge RNAs (sbRNAs) are homologs (or functional analogs) of the Y RNAs. Moreover, analysis of the differential libraries showed that some mature snoRNAs undergo secondary 5' cap modification after processing of the primary transcript, thus suggesting the existence of a wider range of functional RNAs arising from processed and modified fragments of primary transcripts.